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1.  INTRODUCTION 


Calculation  of  the  input  impedance  of  a  dipole  antenna  located 
electrically  close  to  a  plane  dielectric  interface  has  traditionally 
involved  a  Sommerfeld  integral  formulation.  This  is  well  known  to  be 
a  time  consuming  and  difficult  computational  task.  A  good  bibliography 
can  be  found  in  a  book  by  Banos[l].  A  disadvantage  of  the  Scvnerfeld 
formulation  is  that  the  mathematical  manipulations  involved  in  obtain¬ 
ing  a  better  approximation  or  a  better  inteyration  path  seldom  have  a 
clear  physical  explanation.  This  means  that  it  is  often  difficult  to 
decide  what  action  to  take  to  facilitate  the  computation.  Sommerfeld 
solutions  for  finite  length  antennas  are  difficult  to  program  success¬ 
fully  and  the  resulting  program  is  usually  lengthy. 

This  report  presents  a  simple  yet  exact  formulation  for  the  mutual 
impedance  between  two  sinusoidal  dipole  antennas  in  the  presence  of  a 
half  space.  This  formula  can  be  simply  evaluated  numerically  when  the 
dipoles  are  immersed  in  a  lossy  medium,  and  can  be  used  in  a  moment 
method  solution  using  sinusoidal  bases  and  Galerkins  method[2].  The 
resulting  program  is  relatively  compact.  The  mutual  impedance  is  infer¬ 
red  (by  the  Array  Scanning  Method)  from  a  convenient  formula  for  the 
infinite  array  impedance  developed  by  Munk,  et  al [3]  which  contains  a 
sum  of  plane  waves.  The  half  space  is  accounted  for  by  including  a 
plane  wave  reflection  coefficient  for  each  plane  wave.  The  advantages 
of  this  formulation  are  that  it  is  exact,  it  is  simple  to  program,  and 
it  has  a  clear  physical  interpretation. 

The  Array  Scanning  Method  formulation  is  presented  in  Section  II. 
Section  III  discusses  the  result.  The  utility  of  the  method  is  demon¬ 
strated  in  Section  IV  by  the  presentation  of  numerical  results  of  the 
input  impedance  of  dipoles  in  lossy  half  spaces.  Several  different 
cases  are  presented  because  there  are  few  results  of  half  space  calcu¬ 
lations  contained  in  the  literature.  In  Section  V  an  approximation  is 
made  for  low  frequencies  and  conducting  media  which  results  in  a  con¬ 
siderable  saving  in  computer  time.  As  an  example  a  moment  method  solu¬ 
tion  for  the  input  impedance  for  a  dipole  buried  in  the  earth  is  pre¬ 
sented.  The  change  in  the  current  distribution  as  the  dipole  approaches 
to  within  x/150,000  of  the  surface  is  shown. 


II.  THE  MUTUAL  IMPEDANCE  BETWEEN  TWO 
SINUSOIDAL  DIPOLES  IN  TERMS 
OF  PLANE  WAVES 

Munk,  et  al [3]  have  developed  a  convenient  formulation  for  the  free 
space  mutual  impedance  between  a  reference  dipole  located  parallel  to  an 
infinite  two  dimensional  plane  array  of  parallel  sinusoidal  dipoles. 

The  formula  converges  much  faster  on  a  computer  than  does  the  method  in 
which  the  mutual  impedances  of  the  individual  elements  are  sumned.  Re¬ 
ferring  to  Fig.  1,  the  mutual  impedance  between  the  reference  dipole  and 
the  array  when  both  are  in  a  general  medium  of  propagation  constant 
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Fig.  1.  Geometry  of  reference  dipole  and  infinite  array. 


y  =  a  +  jp,  and  all  elements  have  a  sinusoidal  distribution  of  the 
Floquet  type  is 


(1) 


7a  =  2n 

”  dzdx 


P2(e)  e 


JBZso 


where  the  superscript  denotes  an  array  impedance,  and 


is  a  direction  cosine.  The  other  directions  cosines  cx,  c2  are  given 
in  Eqs.  (47)  and  (48)  of  Appendix  A.  The  other  terms  are  1 


y[C0S6£S,-,  -  coshyi]^ 

P2(o)  =  — *-?— % - ?— 

+  y  ]sinh  y*. 
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sin$  + 
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n  =  /\i 71  is  the  characteristic  impedance  of  the  medium,  and  e  and 
<P  are  the  phase  differences  between  the  current  on  adjacent  array  ele¬ 
ments  in  the  z  and  x  directions  respectively.  The  significance  of  the 
summation  indices  will  bs  clarified  in  the  following  section. 

By  applying  simple  Fourier  theory  during  the  derivation  of  Eq.  (1), 
we  instead  obtain  an  expression  for  the  mutual  impedance  between  the 
reference  dipole  and  a  single  array  element  centered  at  the  origin  of 
the  coordinate  system  shown  in  Fig.  1: 


(6) 


fir/8Dz 

•  -n/fiDz 


»/BDx 

Zu  d(sin<j>)d(sine) 

-rr/BDx 


The  details  of  the  derivation  are  given  in  Appendix  A.  The  deri 
vation  of  Eq.  (1)  is  implicit  in  Appendix  A. 


III.  INTERPRETATION  OF  THE  RESULT 

Notice  that  in  Eq.  (6)  that  the  integration  is  done  with  respect  to 
the  variables  sin$  and  sine  which  by  Eqs.  (2),  (47)  and  (48)  control  the 
direction  of  propagation  of  the  plane  waves  emitted  by  the  Infinite  array 
(changing  the  phase  of  the  current  on  the  elements  to  control  the  direc¬ 
tion  of  radiation  is  called  scanning).  Therefore  implicit  in  Eq.  (6) 
is  a  plane  wave  expansion  of  the  spherical  wave  field  emitted  by  a  single 
dipole,  and  the  technique  is  called  the  Array  Scanning  Method  (ASM). 

The  advantage  gained  in  this  approach  is  that  a  clear  physical  inter¬ 
pretation  is  obtained,  viz.  that  the  plane  waves  are  emitted  by  an  in¬ 
finite  array  in  which  the  wanted  antenna  is  a  member  element. 

It  is  instructive  to  examine  the  behaviour  of  the  direction  cosine 
cy.  Figure  2  shows  cy  for  an  array  having  a  spacing  D,  =  3.5A  with  all 
elements  driven  with  equal  phase  (sine  =  sin<j>  =  0),  both  for  free  space. 
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IMAGINARY  AXIS 


REAL  AXIS 


Fig.  2.  Example  of  direction  cosine  cy  for  an  array  in  free  space 
and  in  a  lossy  medium.  For  tne  array,  Dz=3.5*,  sino= 
sinj.=  1^=0.  The  frequency  is  300  MHz. 
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and  for  a  lossy  medium  of  constitutive  parameters  cr=l ,  o  =  .001  mhos/ 
meter.  The  values  are  plotted  for  k0  =  0. 

For  no  loss,  when  nQ  <.  3,  cy  is  real.  This  means  that  the  array 
emits  propagating  waves  in  directions  specified  by  the  direction  cosines 
Cy,  Cy,  cz.  Since  k0  =  0,  cx  =  0  so  that  the  direction  of  radiation  is 
always  in  the  y-z  plane.  The  radiation  associated  with  n0  =  0  occurs  in 
a  direction  specified  by  the  element  phases,  sine  and  sin<t>.  If  sine  = 
s i n<j>  =  0,  the  array  radiates  broadside.  This  is  the  principal  direction 
or  main  lobe.  For  n0  =  1,  2  jr  3,  the  radiation  occurs  in  other  direc¬ 
tions;  these  are  the  grating  lobes. 

When  n0  >  3,  cy  becomes  imaginary,  so  that  the  plane  waves  in  Eq. 

(1)  are  non-propagating  or  evanescent  waves.  The  imaginary  direction 
cosines  specify  the  directions  of  propagation  into  imaginary  space.  The 
evanescent  waves  decay  rapidly  as  the  spacing  d|  from  the  array  increases, 
and  as  n0  increases.  If  is  this  rapid  decay  which  makes  Eq.  (1)  con¬ 
verge  rapidly. 

When  the  medium  is  lossy,  cy  is  complex  and  lies  on  the  curve  shown 
in  Fig.  2.  Thus  the  propagating^ waves  are  attenuated,  and  there  are  no 
purely  evanescent  waves  for  an  infinite  array  in  a  lossy  medium.  How¬ 
ever,  when  nQ  >  3  the  attenuation  factor  becomes  so  large  that  in  prac¬ 
tice  the  waves  decay  in  the  same  way  as  in  the  lossless  case. 

In  Eq.  (6)  it  is  convenient  to  choose  the  array  spacings  Dx  =  Dz  = 
0.5x,  so  that  the  limits  of  integration  become  -1  and  +1.  Figure  3 
depicts  the  region  of  integration  for  the  half-wavelength  spaced  array 
in  free  space.  For  free  space  this  diagram  is  known  as  the  direction 
cosine  plane[4].  It  is  seen  that  Eq.  (6)  corresponds  to  scanning  the 
principal  beam  of  the  array  over  a  cell  consisting  of  all  r^al  space 
(that  within  the  circle)  and  part  of  imaginary  space.  Further,  be¬ 
cause  of  the  form  of  the  summation  in  Eq.  (1)  the  integration  is  effec¬ 
tively  being  carried  out  over  all  space  (in  Fig.  3,  as  the  main  lobe 
scans  the  pictured  square  cell,  the  evanescent  modes  are  scanning  an 
infinite  number  of  identical  square  cells  covering  all  space).  Notice 
that  when  cx  +  cr  >  1,  cy  is  imaginary  and  the  array  does  not  radiate. 
Under  this  condition  the  mutual  impedance  is  purely  reactive.  Notice 
also  that  when  cf;  +  cz  =  1 ,  cy  =  0  and  the  mutual  impedance  is  infinite. 

It  is  this  "infinite  ridge"  on  the  real-imaginary  space  boundary  which 
makes  Eq.  (6)  difficult  to  evaluate  numerically  for  a  free  space  medium. 
However  numerical  experiments  show  that  the  integral  does  converge  when 
the  "pole"  is  approached.  This  verifies  a  conclusion  reached  on  purely 
physical  grounds:  if  the  integral  did  not  converge  the  mutual  impedance 
between  two  dipoles  would  be  infinite  regardless  of  spacing.  However 
when  the  medium  is  lossy,  cy  is  never  zero,  and  it  is  then  straight¬ 
forward  to  evaluate  Eq.  (6rnumerically.  Figure  4  gives  an  c;. ample  of 
the  variation  of  the  mutual  impedance  between  the  reference  element  and 
an  array  in  a  lossy  medium  as  the  scan  angles  0  and  $  are  varied.  The 
resistance  and  reactance  (shown  separately  in  Fig.  4)  are  plotted  for 
the  first  quadrant  of  the  square  cell  depicted  in  Fig.  3.  The  position 
of  the  ridge  is  shown  clearly. 
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IV.  THE  IMPEDANCE  OF  DIPOLES  IN 
A  LOSSY  HALF-SPACE 

When  the  separation  d  of  the  reference  dipole  and  the  array  is  set 
equal  to  the  dipole  radius  a  we  obtain  the  self  impedance[5] .  Next, 
the  mutual  impedance  of  an  image  antenna  in  a  plane  interface  can  be 
computed  simply  by  including  a  plane  wave  reflection  coefficient  for 
each  plane  wave  in  Eq.  (1).  The  input  impedance  of  a  sinusoidal  dipole 
parallel  to  a  plane  interface  is  then  given  by 


(7)  Z,„  -  Z„(a>  *  rn  2H(2d  ) 

0  0 

The  geometry  is  shown  in  Fig.  5.  rn0,k0  is  an  effective  plane  wave 
reflection  coefficient  (i.e.,  including  dipole  polarization  decoupling 
and  surface  waves  when  they  exist)  for  an  array  of  z  directed  dipoles 
parallel  to  a  plane  Interface.  It  is  given  by  [3] 


where  cx>  c  and  c7  are  the  direction  cosines  of  the  incident  plane 
wave,  c  and  are  the  complex  dielectric  constants  of  the  medium 
containing  the  antenna  and  the  half  6pace  respectively,  and  are  given 


by 

(9) 

l  1  -  l  -  ja/uj 

(10) 

^2  = 

8 


MEOIUM  2:«2  ,  <r2 


zla 

MEDIUM  \  <(  ,  <r, 


Fig.  5.  Geometry  of  dipole  parallel  t.o  a  plane  interface. 

The  model  of  Eq.  (7)  obtained  using  a  sinusoidal  dipole  represents  a 
two  segment  moment  method  solution  using  the  sinusoidal  reaction  technique 
[2].  Hence  Eq.  (7)  can  be  used  in  multi  segment  moment  method  solutions 
using  piecewise  sinusoidal  expansion  dipoles  and  Galerkins  method. 

To  demonstrate  the  utility  of  the  ASM  several  examples  of  the  in¬ 
put  impedance  of  a  0.1  m  long,  0.001  m  radius  sinusoidal  dipole  parallel 
to  plane  interfaces  have  been  computed  at  a  frequency  of  300  MHz.  All 
media  are  lossy  dielectrics  (e  >  a)  at  300  MHz.  In  computing  these  re¬ 
sults  the  displacements  X  and  Z  were  zero,  and  advantage  was  taken  of 
the  symmetry  of  Eq.  (6)  about  the  sino  and  s i axes,  i.e.,  integration 
need  only  be  done  for  positive  sine  and  $in$.  The  integration  was  done 
by  the  trapezoidal  rile,  and  the  required  number  of  integration  steps  N 
along  each  axis,  and  the  required  convergence  accuracy  for  the  series 
were  determined  by  numerical  experiment.  This  involved  comparing  the 
mutual  impedance  compjted  using  Eq.  (6)  to  those  computed  by  the  induced 
emf  method  using  closed  form  expressions  for  the  exponential  integrals. 

For  a  frequency  of  300  MHz,  the  following  values  given  in  Table  I 
consistently  gave  accuracies  better  than  2%.  A  reduction  In  the  con¬ 
vergence  accuracy  to  10”^  gave  an  accuracy  of  about  5X. 


TABLE  I 


SUGGESTED  VALUES  FOR  N  AND  CONVERGENCE  ACCURACY 
AT  300  MHz.  VALUES  ARE  FOR  TRAPEZOIDAL  INTEGRATION 


0 

d 

N 

Convergence 

accuracy 

.001 

£  0.2x 

15 

10'4 

.01 

>  0.2x 

15 

lO-J 

£  0.2X 

10 

TO-4 

0.1 

all  d 

10 

o 

■ 

1.0 

all  d 

6 

1 

O 

Since  Eq.  (1)  converges  more  slowly  as  the  distance  d  becomes  small, 
the  computing  efficiency  was  further  improved  by  computing  the  self 
impedance  of  the  dipole  (the  first  term  in  Eq.  (7))  by  the  induced  emf 
method.  The  plane  wave  expansion  was  used  only  to  calculate  the  mutual 
impedance  of  the  image  in  the  interface. 

Figures  6  and  7  show  the  input  impedance  of  the  dipole  when  it 
approaches  an  interface  which  has  a  different  permittivity.  Results 
are  presented  for  the  dipole  on  both  sides  of  the  interface.  Figure  6 
shows  that  when  the  conductivity  a  s  .001  mhos/m  there  is  little  change 
in  the  input  resistance,  but  the  reactance  is  lower  when  the  dipole  is 
in  the  medium  of  higher  permittivity.  This  occurs  because  the  shorter 
wavelength  in  medium  1  means  that  the  dipole  Is  electrically  longer  and 
therefore  is  nearer  to  resonance.  The  smooth  reactance  change  as  the 
dipole  approaches  the  interface  occurs  because  of  the  influence  of  the 
nearby  half-space  on  the  wavelength  as  the  dipole  approaches  the  inter¬ 
face.  When  the  conductivity  of  both  media  is  increased  to  0.01  mhos/m, 
we  observe  in  Fig.  7  that  the  input  resistance  has  increased  significantly. 
This  occurs  because  the  characteristic  impedance  n  for  the  lossier  medium 
has  a  larger  imaginary  part.  Thus  the  evanescent  modes  summed  In  Eq.  (1), 
which  affect  only  the  reactance  when  the  medium  is  lossless  (n  real), 
now  contribute  significantly  to  the  resistance  because  of  the  complex 
n  multiplier  in  Eq.  (1). 

Figures  8-10  show  the  input  impedance  when  the  dipole  approaches 
the  interface  of  a  medium  of  different  conductivity.  When  the  conduc¬ 
tivity  is  low,  as  in  Fig.  8,  there  is  little  change  In  the  wavelength 
and  in  the  real  part  of  n.  so  that  there  Is  little  change  in  the  dipole 
reactance.  The  relatively  large  change  In  the  imaginary  part  of  n  causes 
the  observed  change  in  the  Input  resistance.  In  Fig.  9,  when  the  con¬ 
ductivity  is  Increased  so  that  there  is  a  significant  wavelength  change 
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Input  impedance  of  0.1m  sinusoidal  dipole  in  media  of 
conductivity  0.001,  and  parallel  to  plane  interface 
of  different  permittivity. 
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between  the  media,  a  significant  reactance  change  is  observed.  An  inter¬ 
esting  observation  is  that  practically  all  the  resistance  change  occurs 
when  the  dipole  is  in  the  medium  with  the  lowest  conductivity.  This 
indicates  that  the  medium  of  lowest  characteristic  impedance  absorbs  most 
of  the  power  supplied  to  the  dipole  when  the  dipole  is  close  to  the  inter¬ 
face  but  is  located  in  an  adjacent  medium  of  higher  characteristic  imped¬ 
ance.  This  effect  is  even  more  pronounced  in  Fig.  10  where  the  two  media 
differ  in  conductivity  by  two  orders  of  magnitude. 

Because  of  their  simplicity  dipole  antennas  make  good  subterranean 
radar  antennas[6].  Consequently  the  characteristics  of  dipoles  buried 
in  the  ground  near  the  earth-air  interface  are  of  interest.  Figure  11 
shows  how  the  input  impedance  of  a  0.001  m  radius,  0.1  m  long  dipole 
varies  with  different  burial  depths  in  earth  having  constitutive  para¬ 
meters  t'r  =  4,  a  =  .01  mhos/m.  Three  curves  are  shown.  One  is  a  direct 
application  of  Eq.  (7).  For  comparison  results  obtained  from  a  four 
segment  moment  method  using  sinusoidal  bases  and  Galerkins  method  with 
Eq.  (7)  to  obtain  the  impedance  matrix  elements  are  also  presented. 

Note  that  for  the  electrically  short  dipole  considered  good  accuracy  is 
obtained  with  a  two  segment  solution.  This  is  further  confirmed  by  the 
eight  segment  solution  by  Richmond[7]  for  a  homogeneous  medium  (i.e., 
d  =  “). 

Also  presented  in  Fig.  11  is  a  comparison  with  a  Sommerfeld  inte¬ 
gral  solution  by  program  WFLLL2B[8] .  It  is  emphasized  that  with  this 
program  21  segments  were  required  to  obtain  the  results  presented,  which 
is  equivalent  of  105  segments  per  wavelength  in  the  earth.  In  comparison 
it  is  seen  thal  good  results  are  obtained  for  this  antenna  (which  is  0.2a 
long  in  the  earth)  by  using  only  two  segments  and  the  ASM.  It  is  noted 
also  that  with  WFLLL2B  and  21  segments  there  is  still  some  difficulty  in 
matching  the  result  for  a  homogeneous  medium  with  those  for  the  interface. 
Although  the  reactances  obtained  by  the  two  methods  differ  slightly,  the 
results  are  consistent. 

The  behaviour  of  the  impedance  in  Fig.  11  is  consistent  with  that 
observed  for  a  dipole  approaching  a  medium  of  lower  conductivity  shown 
in  Figs.  9  and  10.  Note  that  results  could  not  be  computed  for  the 
dipole  in  air  with  the  interface  present.  A  moment  method  solution[7] 
using  8  segments  in  free  space  yielded  a  value  of  1.69  -  j 1 036  for  the 
input  impedance.  Note  that  the  effect  of  conductive  earth  contact  is 
to  greatly  increase  the  input  resistance.  As  noted  before,  the  react¬ 
ance  change  is  primarily  due  to  the  change  in  the  electrical  length  of 
the  antenna  as  the  interface  is  approached.  The  resistance  is  largely 
unaffected  because  of  the  power  input  to  the  antenna  being  absorbed  by 
the  earth. 

The  computation  time  is  observed  to  vary  approximately  as  d"^.  Con¬ 
sequently  the  results  for  small  d  are  relatively  time  consuming.  Table 
II  shows  the  computer  time  required  to  calculate  the  results  shown  in 
Fig.  6  for  two  convergence  accuracies.  The  computer  used  was  a  Datacraft 
6024,  but  the  times  have  been  presented  for  an  IBM  370/168  computer  by 
running  one  example  and  converting  all  times  by  the  speed  ratio  thus 
obtained. 
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TABLE  II 


IBM  370/168  COMPUTER  TIME  REQUIRED  TO  COMPUTE 
THE  RESULTS  PRESENTED  IN  Fig.  6 


distance  from 
interface  d 

Seconds  of  CPU 
time  for 

convergence  =  10-3 

Seconds  of  CPU 
time  for 

convergence  =  10-/* 

.5  m 

in 

CO 

» 

0.87 

.25  m 

1  .05 

1  .23 

.1  m 

3.37 

4.27 

.05  m 

8.41 

12.4 

.01  m 

28.6 

46.3 

.005  m 

45.6 

103 

.001  m 

117 

477 

For  d  =  .001  m,  the  impedance  calculated  with  a  convergence  accur¬ 
acy  of  10-3  was  4,697  -  j!48.4r..  When  the  convergence  accuracy  was  in¬ 
creased  to  TO-4  the  calculated  impedance  was  4.285  -  jl43.8ii  It  is 
doubtful  if  the  improvement  in  accuracy  justifies  the  large  increase  in 
computer  time.  Large  amounts  of  computer  time  are  required  when  d  is 
less  than  about  x/1000. 

In  Table  III  a  CPU  time  comparison  is  made  between  the  ASM  and 
the  Sommerfeld  program  WFLLL2B[8].  Since  the  ASM  is  similar  to  a  Somner- 
feld  approach[3]  no  dramatic  differences  in  CPU  time  are  observed.  It 
is  apparent  however  that  far  fewer  segments  are  needed  with  the  ASM 
because  of  the  superiority  of  the  piecewise  sinusoidal  current  expansion 
for  wire  antennas.  It  should  also  be  noted  that  the  pilot  subroutine 
used  trapezoidal  integration.  A  significant  improvement  is  likely  if  a 
more  efficient  integration  procedure  is  used. 

It  is  also  worth  pointing  out  that  our  pilot  subroutine  which  com¬ 
putes  the  values  of  the  integrand  contains  only  65  FORTRAN  statements. 

A  further  13  statements  are  added  to  calculate  r„  l  ,  and  the  integra¬ 
tion  subroutine  contains  70  statements. 
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TABLE  III 


IBM  370/168  CPU  TIME  IN  SECONDS  REQUIRtD  TO 
COMPUTE  THE  RESULTS  PRESENTED  IN  Fig.  11 


d 

(meters) 

Array  scanning  method' 

Program 

WF-LLL-2B 

21  segments 

4  segments 

2  segments 

.25  m 

5.6 

1 .2 

68 

1  m 

17.5 

4.2 

59 

.05  m 

56.3 

12.3 

79 

.01  m 

341 

45 

105 

.005  m 

526 

101 

136 

V.  THE  IMPEDANCE  OF  DIPOLES  IN  A 

CONDUCTING  HALF-SPACE 

Subterranean  electromagnetic  pulse  radar  systems  for  probing  many 
hundreds  of  meters  into  the  earth  must  necessarily  operate  at  low  fre¬ 
quencies^].  Physically  long  antennas  lying  on  the  surface  or  buried 
only  a  few  centimeters  from  the  surface  are  required.  This  close 
electrical  spacing  to  the  surface  threatens  the  viability  of  this  method 
for  calculating  the  input  impedance  cf  these  antennas.  Fortunately 
advantage  can  be  taken  of  the  low  frequency  characteristics  of  the  earth. 

At  low  frequencies  where  the  earth  is  classed  as  a  conductor  (a  = 
tf),  xhe  higher  attenuation  constant  assists  convergence  (rate  of  con¬ 
vergence  for  large  n0  and  k0  in  Eq.  (1)  depends  upon  the  real  part  of 
yc„,  which  varies  in  a  rather  complicated  way  with  frequency  and  con¬ 
stitutive  parameters).  Further  economies  in  computer  time  are  possible 
at  low  frequencies  by  reducing  the  number  of  integration  steps  N.  This 
is  due  to  the  increased  attenuation  between  the  array  elements  as  the 
wavelength  becomes  large  (which  means  that  the  array  spacing,  which  we 
set  at  0.5x  in  our  numerical  calculations,  becomes  large). 

In  this  section  calculations  of  the  input  impedance  of  a  1000  m 
dipole  in  earth  of  constitutive  parameters  er  =  4,  o  =  .01  mhos/meter 
are  presented  for  a  frequency  of  100  Hz.  It  was  possible  to  calculate 
results  with  the  antenna  as  close  as  2  cm  to  the  interface  with  reason¬ 
able  expenditure  of  computer  time  (2  cm  is  about  1/150,000  of  a  wave¬ 
length  in  the  earth).  For  this  example  good  results  were  obtained  when 
the  dipole  was  within  O.lx  (^  300m)  of  the  interface  by  evaluating 
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only  the  integrand  of  Eq.  (6)  (i.e.,  N  =  0)  with  sino  =  sin<j,  =  0.  This 
means  that  the  mutual  impedance  of  the  array  is  almost  identical  to  the 
mutual  impedance  of  a  single  element.  This  occurs  because  at  100  Hz 
the  array  is  spaced  x/2  r  1581  meters  between  elements,  and  the  high 
loss  means  that  only  the  fields  of  the  center  element  are  significant 
providing  the  reference  dipole  is  located  near  the  center  element. 

Figure  12  plots  the  error  in  the  mutual  impedance  of  two  parallel 
1000m  sinusoidal  dipoles  in  an  infinite  medium  of  i.r  =  4,  o  =  .01  mhos/ 
meter  when  calculated  using  this  approximation.  The  induced  emf  method 
was  used  as  a  reference.  The  error  is  seen  to  increase  as  the  separation 
of  the  dipoles  is  reduced.  The  error  also  increases  when  the  separation 
increases  beyond  300m:  this  occurs  because  the  spacing  between  the  refer¬ 
ence  dipole  and  the  adjacent-to-center  elements  of  the  array  is  less 
than  an  order  of  magnitude  smaller  than  the  spacing  of  the  reference 
dipole  and  the  center  element.  This  error  can  be  reduced  by  performing 
the  integration  in  Eq.  (6),  i.e.,  N  >  0.  However  the  integration  was 
observed  to  have  little  effect  on  the  error  at  close  spacings,  supporting 
the  assumption  that  the  mutual  impedance  is  dominated  by  the  fields  of 
the  center  element  of  the  array.  The  computer  time  required  is  indicated 
on  Fig.  12:  the  saving  is  seen  to  be  roughly  50x  by  comparison  with 
Table  II  for  similar  electrical  spacings. 

Using  the  array  approximation,  the  input  impedance  of  the  1000m 
dipole  buried  in  the  earth  was  calculated.  The  results  are  shown  in 
Fig.  13.  The  convergence  accuracy  chosen  was  0.5  x  10“3.  Since  the 
induced  emf  method,  which  is  used  to  calculate  the  self  impedance,  is 
assumed  accurate,  and  Fig.  1?  gives  the  expected  accuracy  of  the  mutual 
impedance  of  the  image  in  the  interface  assuming  perfect  reflection,  the 
overall  accuracy  of  the  calculations  presented  in  Fig.  13  can  be  estimated. 
The  error  in  the  reactance  is  estimated  to  be  less  that  1%  for  all  d,  and 
the  error  in  the  input  resistance  is  estimated  to  be  less  than  2%  when 
d  =  1m. 


When  d  =  0.1m,  the  error  in  the  resistance  exceeds  8%.  It  is  inter¬ 
esting  to  note  that  when  d  <  0.1m  the  error  in  the  mutual  impedance  com¬ 
putations  acts  to  make  the  mutual  (and  hence  self)  impedance  constant. 

This  effect  was  also  noticed  when  preparing  Fig.  12,  hence  it  is  not  a 
property  of  the  half-space  (i.e.,  although  the  error  increases  as  the  spac¬ 
ing  decreases,  the  mutual  impedance  calculated  by  the  ASM  is  observed  to 
become  independent  of  spacing). 

The  economies  in  computer  time  (indicated  by  the  CPU  times  included 
in  Fig.  13)  obtained  by  the  array  approximation  described  in  this  section 
make  this  technique  useful  for  more  accurate  analysis  of  dipoles  by  the 
moment  method.  As  an  example,  the  input  impedance  of  the  1000m  dipole 
obtained  from  a  four  segment  (equal  lengths)  moment  method  solution  using 
sinusoid;^  bases  and  Galerkins  method  is  presented  in  Fig.  14.  Note  the 
improvement  in  accuracy  for  small  d,  and  the  appearance  of  the  "kink"  near 
d  =  c°,  commonly  observed  in  the  lossy  dielectric  results.  Results  for  a 
higher  series  convergence  accuracy  (10-^)  are  given  in  Fig.  15.  As  concluded 
in  the  previous  section,  the  improvement  in  accuracy  does  not  warrant  the 
large  increase  in  computer  time. 
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Fig.  13.  Input  impedance  of  1000m  long  sinusoidal  dipole 
buried  in  a  conducting  half  space  (earth).  The 
computation  times  in  seconds  for  an  IBM  370/168 
computer  are  given  in  parenthesis. 


It  is  interesting  to  note  from  Fig.  15  that  when  the  dipole  approaches 
the  air  interface  from  the  conducting  medium  that  the  only  significant 
change  occurs  in  the  input  resistance.  The  results  presented  indicate 
that  if  the  dipole  was  brought  from  d  =  »  into  close  contact  with  the 
interface  then  the  input  resistance  would  approximately  double.  This 
confirms  the  simple  physical  argument:  there  is  only  one  conducting  half 
space  to  shunt  the  antenna  so  that  the  input  resistance  is  doubled. 
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Fig.  14.  Four  segment  sinusoidal  Galerkin  moment  method  solution 
for  input  impedance  of  1000m  long  dipole  burled  in  a 
conducting  half-space  (earth).  The  series  convergence 
accuracy  was  0.6  x  10'3.  IBM  370/168  computation  times 
in  seconds  are  given  in  parenthesis. 
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Fig.  15.  Four  segment  sinusoidal  Galerkin  moment  method  solution 
for  input  impedance  of  1000m  long  dipole  buried  in  a 
conducting  half  space  (earth).  The  series  convergence 
accuracy  was  10  .  IBM  370/168  computation  times  in 

seconds  arc-  given  in  parenthesis. 
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Mul tisegment  moment  method  solutions  with  small  segments  are  more 
accurate  than  simple  two  segment  solutions  because  the  current  distri¬ 
bution  is  more  accurately  represented.  With  a  two  segment  solution  it 
is  assumed  that  the  current  distribution  remains  unchanged  as  the  dipole 
approaches  the  interface.  The  accuracy  of  this  assumption  can  be  gauged 
from  the  four  segment  solution.  Table  IV  shows  the  current  relative  to 
the  terminal  current  at  a  point  halfway  between  the  center  terminals  and 
the  dipole  tip  as  the  dipole  aporoaches  the  interface.  If  the  medium 
containing  the  dipole  was  lossless,  the  magnitude  of  the  current  would 
have  been  slightly  greater  than  0,5  (since  the  sinusoidal  shape  of  the 
distribution  begins  to  show  when  the  antenna  is  0.3x  long)  and  the  phase 
near  zero.  The  attenuation  of  the  medium  acts  to  reduce  the  magnitude 
of  the  current  along  the  antenna,  and  causes  the  negative  phase  angle. 

The  physical  explanation  is  that  the  reflection  of  the  input  wave  from 
the  tips  of  the  antenna  is  attenuated,  and  when  summed  with  the  outgoing 
wave  results  in  a  current  amplitude  with  a  negative  phase  angle.  Notice 
in  Table  IV  that  as  the  antenna  comes  close  to  the  surface  the  immediate 
environment  of  the  antenna  becomes  less  lossy,  thereby  reducing  the  phase 
and  increasing  the  magnitude  of  the  current  at  the  midpoint  of  each  side 
of  the  antenna . 


TABLE  IV 

AMPLITUDE  OF  THE  CURRENT  RELATIVE  TO  THE  TERMINAL 
CURRENT  250m  FROM  THE  TERMINALS  OF  A  CENTER  FED 
1000m  DIPOLE  IN  CONDUCTING  EARTH  PARALLEL 
TO  AN  AIR  INTERFACE 


distance  from 
interface  d 
(metres ) 

Relative  amplitude 
of  current 

300 

0 .47 0/ -15  .7° 

100 

0 .46U>-15-.£u 

60 

0.463/ -15 .5° 

40 

0. 464/ -J  5^1° 

20 

0.464/ -14  .6° 

10 

0 .465/jJL2ju£° 

5 

0.467/JXJ° 

2 

0.470 L12^A° 

1 

0.471/^1^2° 

0.5 

0.473/-) 1  .4° 

0.1 

0 .476/-10  .6° 

0.05 

0 .476/  -10  .4° 

0.02 

0.477/-10.30 
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VI.  CONCLUSIONS 


A  simple  alternative  formulation  of  the  classical  Sommerfeld  pro¬ 
blem  called  the  Array  Scanning  Method  (ASM)  has  been  applied  to  linear 
sinusoidal  antennas  in  a  lossy  half  space.  The  advantage  of  the  method 

are : 

1.  It  is  simple  to  program  and  the  resulting  program  occupies 
significantly  less  storage  than  does  a  typical  Soimerfeld 
program. 

2.  The  formulation  has  a  clear  physical  interpretation  which 
greatly  facilitates  its  use. 

It  is  shown  that  at  low  frequencies,  the  input  resistance  of  a 
dipole  buried  in  the  earth  approximately  doubles  when  it  is  moved  from 
large  depths  to  a  position  near  the  interface.  This  confirms  the  simple 
physical  argument:  there  is  only  one  conducting  half-space  to  shunt  the 
antenna  so  that  the  resistance  is  doubled. 

Although  the  numerical  evaluation  of  the  integral  is  not  as 
straightforward  when  the  antenna  is  in  free  space,  the  method  has  been 
shown  to  converge  in  this  case,  and  with  a  suitable  integration  pro¬ 
cedure  could  be  applied  to  the  analysis  of  linear  antennas  above  the 
earth. 
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APPENDIX  A 


MUTUAL  IMPEDANCE  OF  PARALLEL  SINUSOIDAL  DIPOLES 
BY  THE  ARRAY  SCANNING  MFTHOD 


Consider  a  reference  dipole  located  parallel  to  an  infinite  plane 
array  in  a  lossy  medium  with  c  and  o,  as  shown  in  Fig.  1.  The  displace¬ 
ments  from  the  x  axis  are  X  and  Z.  Following  Munk,  et  al[3]  the  current 
on  the  element  in  the  kth  column  and  nth  row  of  *he  array  is  given  by 

j6nDxsin6  jrtkD,sin<t> 

<">  *kn  '  1  e 


where  I  =  sinhy(8.-|z  |  )/sinhy£  is  the  current  on  the  centre  element, 
y  =  a+ j 6  is  the  propagation  constant,  and  e  and  <t>  are  phase  differences 
in  the  current  on  adjacent  elements  in  the  2  direction  and  x  direction 
respectively. 

The  mutual  impedance  between  the  reference  dipole  and  the  array  Zu  can 
be  obtained  by  summing  the  mutual  impedances  of  each  element  of  tne  array 

”  “  j6nDxsine  jekD7sin<t> 

on  zj  ■  l  l  z0>k„ « 

k=-»  n--» 

1Q  kn  is  the  mutual  impedance  between  the  reference  dipole  and  the  knth 
element  of  the  array,  and  is  given  by[5] 


(13) 


2o,kn 


-.1 

Airsinh^yt 


e“Yri 


+ 


e~yr2 

r2 


-  2  coshyH 


sinhy( l- | 2 | )dz 


where  n  is  the  characteristic  impedance  of  the  medium,  and 

(14)  r2  =  (kOx  -  X)2  +  d2  +  (nDz  -  Z  -  z)2 

(15)  r2  =  (kDx  -  X)2  +  d2  +  (nDz  -  Z  -  z  +  i)2 

(16)  r2  =  (kDx  -  X)2  +  d^  +  (nDz  -  Z  -  z  -*)2 

as  shown  ir.  Fig.  16. 
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Fig.  16.  The  geometry  involving  the  reference  dipole  at 
(X,-d,Z)  and  the  kn'th  dipole  of  the  array. 


Equation  (12)  can  be  written  as 


07)  Z*  =  l  he 

k=-c° 


jiBkOzSin* 


where 


(is)  z,  =  — y~ 

4-nsinh  yi. 


”  f*-  p'Yrl  e' 

y  § L  +  *L 

L  «  r, 

n=-°°  -t  L  I 


2e'Yro 


COShyt 


jBnD-sino 


s i  nh'ft  n-  |  z  |  )d z 


Reversing  the  order  if  integration  and  surmation  in  Eq.  (18)  gives 
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(19) 


Zi  = 


-n 


4irsinh2Yt  i 


£  00 

[  sinhy (a- 1 z  | )  l 
1-1  n=-» 


-yr-j  -yrD  -tr„ 

1  +  e _ £  2  e _ 0 


rl 


r2 


jenDzsine 


dz 


By  combining  the  Poisson  sum  formula 

jnwt 


(20) 


Y  F(nu)e 

n=  -00 


T  l  f(t+n0T)  ; 


no_'“ 


and  the  frequency  shift  theorem 


Fourier 

(21)  F(w-W)  «• - >  f(t)  e 


jWt 


T  = 


2v 


where  W  is  a  radian  frequency  shift, we  obtain 

00  .  ®  jW(t+n  T) 

(22)  l  F(«-W)  eJnu)t  =  T  T  e 

n=-®  no*-“ 


f(t+n0T) 


The  summation  in  Eq .  (18)  can  be  done  by  applying  Eq.  (22)  with 
the  Fourier  Transform  pair 


-Yja^U-oq)2  Four1er  j  t  T----T 

(23)  g - --  ? - >  e  '  ’  4  H(2)  (aj-y2-t2) 


>  a2+(u-i 


«1 ) 


and  by  making  the  following  substitutions 


(24) 

(25) 

(26) 


OJ  =  D2 

t  =  e  sine 

z-i  for  r-j 
z  for  rn 


J1 


z+t  for  r, 


coshvt 
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(27) 

(28) 


U  »  Z 


a2  =  (kDx  -  X)2  +  d 

(29)  T  =  ~ 

1Jz 

Equation  (18)  then  becomes 

(30)  Z1  = 


2Dzsi 


:i— -  [  sinhY(t-|z|)  l 

inhS*  n0= 


where  se  is  defined  in  Eq.  (5). 

After  doing  the  integration  we  get 

jsZs 


l)  •  ^P-  l  P2(e)e 

U  t  n  —  —  rr, 


z  No¬ 


where  P^(e)  is  defined  in  Eq.  (3). 

Zj  is  the  mutual  impedance  between  the  reference  dipole  and  an  in¬ 
finite  column  of  dipoles  spaced  Dz  metres  apart  in  the  z  direction.  The 
mutual  impedance  between  the  centre  dipole  only  and  the  reference  dipole 
can  now  be  inferred  from  the  column  array  impedance  Eq.  (31)  by  using 
the  Fourier  series  for  a  periodic  function: 
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Making  the  substitutions  given  in  Eqs.  (24)  -  (29),  we  obtain 


.tt/ SD, 

(35)  ZQ0  =  d(sine) 

;  -7f/6Dz 


ZQ  0  is  the  mutual  impedance  between  the  reference  dipole,  and  the 
centre  element  of  the  z  oriented  column  of  dipoles  located  at  z=0.  Con¬ 
sider  Eq.  (35)  when  the  element  spacing  02  =  x/2.  The  integration  limits 
then  correspond  to  -v/2  <_  e  <_  n/2.  Physically  Eq.  (35)  moans  that  the 
phase  of  the  current  on  the  elements  adjacent  to  the  centre  element  is 
varied  by  tt  radians,  the  next  elements  by  2*  radians,  and  so  on,  so  that 
when  the  integration  is  done  their  mean  field  is  zero  at  the  reference 
dipole.  Thus  all  elements  but  the  centre  element  are  "phased  out".  A 
time-frequency  analogy  is  that  the  d.c.  term  is  the  average  value  of  the 
time  signal . 

By  substituting  Eq.  (35)  into  the  expression  for  the  mutual  imped¬ 
ance  of  the  infinite  array  Eq.  (12)  we  obtain  the  mutual  impedance  of 
the  reference  dipole  and  an  infinite  number  of  dipoles  spaced  along  the 
x  axis  Z0>k0 


(36) 


2  o.ko 


•n/eD2 

■  -n/eDz 


P2(e) 


jeZs 

e 


6 


(*j 


?  o  o. 


d(sinu ) 


jekD2sintt) 

e 


Reverse  the  order  of  summation  and  integration  to  obtain 


(37) 


-o.ko 


(■t/BD, 


i  ej8ZS°  l 

=  -■»  k=-«- 


-ti/BDz  n0=-“ 


jekDzs  i  n4> 
e  d(sine) 


The  k  summation  can  be  summed  by  again  using  Eq.  (22)  with  the 
Fourier  transform  pair 


(38) 


-jd|[-Y2-B2s2]-t2 

ij[-r2-82s§3-t2 


Fourier  H^2^( 


[-y2-B2s2]  JuAd2'  ) 


The  necessary  substitutions  are 


! 

h 


(39) 

T  = 

Dx 

(40) 

t  =  hsin? 

(41) 

““  Dx 

(42) 

W  =  X, 

so  that 

(43) 

,  .  r/sDz 

0'k0"DxDZ  i-./SD 

2  0 


00  ^  JhZS  " 

l  PL(f)e  0  l 

*o=-° 


jeXsj 


-jdJ-Y2“;,2(s|+s^} 

J  -y2-f2(s?+s*) 

V  V 


.  d(sine) 


where  s  is  given  in  Eq.  (5). 

<p 

By  again  applying  Eq.  (34)  with  Eqs.  (39)  -  (42)  we  eliminate  all  but 
the  centre  element.  The  result  is 


d(sin^)a(sine) 


(44) 

-ti/bDz  ,  TI  /  6  Dx 

0,00  -7T/BDZ  J  -TT/6D, 

where 

(45) 

2«  *  6J|n  l  pZ(e)e 

jeZs  a 


l  e‘ 


JBXS 


4>  e 


-Jdj -y2_32(s2^s2) 


ko’'° 


Equation  (45)  is  the  mutual  impedance  between  the  reference  dipole  and 
an  infinite  array,  and  Eq.  (44)  shows  how  the  mutual  impedance  of  a 
single  element  can  be  inferred  from  the  array  impedance  by  phasing  out 
all  elements  but  the  centre  element. 
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1 


The  important 
waves.  Recognizing 


feature  of  Eq.  (45)  is  that 
this  it  car  be  written  as 


it  contains  a  sum  of  plane 


(46) 


7  _  _2n_ 

ZH  '  0  0 
x  z 


P2(o)e 


e 


JBxs 

e 


-vdcv 
e  J 


where  the  direction  cosine  Cy  for  the  plane  waves  emitted  by  the  array 
is  defined  in  Eq.  (2) . 

To  obtain  Eq.  (46)  it  was  necessary  to  use  =  -j  so  that  the 
waves  emitted  by  the  array  attenuate  as  they  propagate  away  from  the 
array.  In  Eq.  (2)  the  required  root  is  the  principal  root.  The  other 
direction  cosines  are 


(47) 


cx  = 


y 


(48) 


J6S„ 
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